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ABSTRACT 


/? 


The  rational  modeling  and  empirical  correlations  used  to  build 
a  comprehensive  computer  code  for  simulating  general  aix-degree-of- 
freedom  motions  of  missile  debris  fragments  are  described.  The 
approach  ia  deterministic  in  that  a  number  of  possible  generic 
fragment  shapes  were  defined,  methods  were  selected  to  describe 
the  aerodynamic  loads  on  these  shapes,  and  the  results  were 
incorporated  In  a  six-degree-of-freedom  trajectory  program.  The 
method  chosen  is  simple  enough  to  avoid  large  computation  times 
and  yet  represents  coning  and  tumbling  conditions  as  well  as 
trimmed  flight. 


> 


INTRODUCTION 


NAVSEA  has  an  extensive  program  concerned  with  point  defense  of  targets 
against  Incoming  missiles.  When  a  defensive  weapon  system  damages  an  attacking 
missile,  the  trajectories  of  the  debris  fragments  remain  of  intarast,  since 
they  may  hit  the  target  or  nearby  areas  and  cause  significant  damage.  Clearly, 
the  possibility  of  such  an  occurrence  increases  with  decreasing  intercept 
distance  and  with  increasing  attacker  velocity.  To  determine  minimum  intercept 
distances  for  a  given  level  of  probable  damage,  it  is  necessary  to  estimate  the 
trajectories  of  the  various  portions  of  debris  following  a  missile  breakup.  \ 
This  requires  calculation  of  the  aerodynamic  forces  and  momenta  over  0-180° 
angle  of  attack  range  and  a  wide  range  of  flight  speeds,  angular  rates  and 
acceleration  conditions.  The  methods,  of  course,  are  also  applicable  to  range 
safety  calculations. 

Under  contract  with  NAVSEA,  Nielsen  Engineering  &  Research,  Inc.  (NEAR) 
has  developed  engineering  methods  for  prediction  of  the  aerodynamic  character¬ 
istics  of  missile  debris  fragments  in  slx-degree-of-f readom  motion.  The 
purpose  of  this  paper  is  to  outline  those  methods,  to  describe  the  resulting 
computer  program,  DEBRIS,  and  to  give  representative  sample  calculations. 


OVERALL  APPROACH 


To  be  useful,  the  methodology  adopted  for  the  simulation  of  six-degree-of- 
freedom  trajectories  of  missll  ’  agments  must  encompass  realistic  ranges  of 
fragment  shapes,  attitudes  am  .  itiee.  Although  a  very  wide  range  of 
fragment  shapes  is  possible,  this  work  is  concerned  with  a  set  of  generic 
shapes  which,  depending  on  conditions,  could  be  stable,  tumble  or  trim  to  a 
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nonzero  angle  of  attack.  Since  the  various  shapes  can  tumble  or  trim  depending 
on  the  relative  locations  of  the  fragment  center  of  gravity  and  aerodynamic 
center  of  pressure,  the  methodology  has  been  designed  to  compute  the  fragment 
aerodynamic  characteristics  for  pitch  and  yaw  angles  up  to  180®  and  roll  angles 
up  to  360*.  Any  spaed  is  allowed  except  that  the  empirical  input  for  the 
methods  described  was  developed  for  the  Mach  range  0.8  to  3.0. 

the  fragment  shapes  fox  which  tha  method*  may  ba  used  are  shown  in  Figure  1. 
Note  that  there  ere  three  beaic  shapes!  (1)  a  cylindrical  body  alone  with 
circular  cross  section  end  with  or  without  pointed  notaj  (2)  body  with  one  aet 
of  identical  fine;  (3)  body  with  two  sate  of  fine,  each  eat  having  Identical 
fins.  Each  finned  section  may  have  up  to  four  fins  set  at  arbitrary  circum¬ 
ferential  positions  on  tha  body. 

Except  for  aerodynamic*  methodology,  the  beaic  framework  for  the  computer 
program  waa  provided  by  the  Science  Applications,  Inc.  (SAI)  coda  known  aa  MAT6. 
The  SAI  program  consists  of  a  main  program,  a  Rungs-Kutta  integration  subroutine, 
a  routine  to  perform  coordinat*  t rant format ion*  and  calculate  the  derivatives 
of  tha  aquations  of  motion,  a  table  look-up  subroutine  for  computing  aarodynamio 
forces  and  momenta  and  an  autopilot  algorithm.  A  description  of  the  methodology 
and  a  list  of  tha  program  is  given  in  Reference  1.  For  tha  purposes  of  this 
work  the  aubroutin*  supplied  by  SAI  for  tha  computation  of  aerodynamic  forces 
and  moments  has  bean  replaced  by  e  new  routine,  GENERIC,  and  it*  satellite 
routines.  The  new  program  ha*  been  named  DEBRIS. 

Subroutine  GENERIC  and  it*  satellite  routines  rapreeent  the  methodology 
derived  during  the  present  work.  They  compute  the  aerodynamic  forces  end 
moments  acting  on  a  debris  fragment  given  tha  kinematic  state  of  the  fragraont 
and  tha  local  propartlaa  of  tha  atmosphere.  Mott  of  the  computations  era  done 
by  satellite  routines  with  GENERIC  acting,  in  affect,  aa  a  driver  program  for 
the  force  and  moment  calculation*.  For  a  typical  case,  GENERIC  perform*  the 
following  steps  aa  it  builds  up  the  loade  for  e  complete  fragment! 

(e)  Compute  forces  end  momenta  acting  on  the  fuselage  or  body  of  the 
debrie  fragment  as  if  no  fins  were  present. 

(b)  Compute  two  velocity  components  at  the  area  centroid  of  each  fin  In 
the  wing  section.  The  firet  component  is  parallel  to  the  body  axis. 

The  second  la  perpendicular  to  the  plane  the  fin  would  occupy  if  It 
were  not  deflected.  These  velocity  components  include  the  Increments 
due  to  rigid  body  rotation  about  the  fragment  center  of  gravity. 

(c)  Compute  the  "equivalent"  angle  of  attack  for  each  fin  in  the  wing 
section. 

(d)  Compute  normal  force  and  canter  of  pressure  for  each  wing  fin.  The 
equivalent  angles  of  attack  computed  above  era  uaad  to  determine 
tha  fin  normal  forces  baaed  on  a  wing-alone  correlation.  A  similar 
approach  Is  used  to  obtain  the  fin  centara  of  prtesurt. 

(e)  Compute  overall  forces  end  moments  due  to  wing  section.  The  methods 
of  Pitts,  Nielsen,  end  Kaattarl^  are  used  to  determine  body  carryover 
loads  due  to  the  presence  of  the  wing  section. 
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(f)  Compute  wing-tall  interference.  At  thla  point*  an  estimate  is  made  of 
the  interference  effects  of  the  wing  section  on  the  tail  section. 

This  is  done  by  estimating  an  equivalent  angle  of  attack,  for  each 
tail  fin  due  to  the  trailing  vortices  shed  by  the  wing  section. 

First*  an  estimate  of  the  strength  and  location  of  those  vortices 
over  Che  tall  section  is  made.  Than  the  methods  of  Reference  2  are 
used  to  get  the  overall  dovnwash.  Finally*  the  appropriate  compo¬ 
nents  of  that  downvash  for  each  fin  are  computed. 

(g)  Steps  (b)  through  (e)  are  repeated  for  the  tall  section. 

The  program  is  also  capable  of  estimating  thrust  effects*  eide  forces  and 
moments  due  to  a  transverse  jet,  and  out-of-plane  aide  forces  and  momenta  due 
to  asymmetric  vortex  shedding.  To  gain  confidence  in  the  computations  for 
debris  fragments,  the  user  can  compute  Che  static  forces  and  moments  for  the 
complete  vehicle  configuration  (before  break-up)  and  compare  those  reaults 
with  available  data.  Complete  details  of  the  methodology  are  available  in 
Reference  3. 

In  the  rest  of  the  paper*  the  procedures  used  for  computing  body-alone 
tore' i  and  moments,  fin  equivalent  angle  of  attack,  fin  forces  and  wing-tall 
Interference  are  given.  The  paper  is  concluded  with  example  calculations  and 
a  summary. 


PROCEDURES  FOR  BODY-ALONE  METHOD 

The  forces  and  moments  acting  on  a  slender  body  without  fins  are  obtained 
from  Blender  body  theory  augmented  by  the  crossflow  theory  of  Allen  (Reference 
4,  Ch.  4).  The  effects  of  rigid  body  rotation  about  tha  body  canter  of  gravity 
have  been  Included.  However,  acceleration  terms  (l.s.  terms  dependent  on  u,  v, 
w,  p,  $,  f)  have  been  neglected.  Two  physical  affects  are  modeled) 

(1)  loading  due  to  instantaneous  values  of  the  crossflow  velocitieij 

(2)  loading  due  to  apparent  mass  effects  which  arise  when  there  le  a 
variation  of  crossflow  velocity  with  axial  position  on  the  body 
axis  (Reference  4,  Ch.  10). 

Whan  crossflow  velocities  are  high  enough  compared  to  tha  axial  velocity, 
flow  separation  on  the  leeward  side  of  the  body  will  occur.  The  effects  of 
separation  on  the  loading  of  type  (1)  above  are  accounted  for  by  including 
crossflow  drag.  It  is  assumed  that  the  type  (2)  loading  above  is  unaffected 
by  separation.  This  is  equivalent  to  assuming  that  the  apparent  mass  for 
viscous  flow  over  a  cross  section  of  the  body  is  the  seme  as  for  lnvlscld  flow. 
Since  the  apparent  mass  is  primarily  dependent  on  the  body  cross  sectional 
dimension  normal  to  the  crossflow  and  since  the  separated  flow  affects  this 
dimension  only  slightly,  the  assumption  appears  to  be  reasonable. 

The  appropriate  expressions  for  the  forces  and  moments  have  been  developed 
previously  by  Goodwin  at  al.  (Reference  3,  eqs.  (46)- (49)  and  (31)- (34))  and 
are  presented  below. 
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The  coordinate  system  used  Is  ohown  in  Figure  2.  It  has  been  assumed  that 
is  independent  of  position  along  the  body  axia.  This  assumption  requires 
that  the  variations  along  the  body  of  crossflow  Mach  number  and  Reynolds 
number  not  be  significant  or 

r (x  -  xa>  <<  v 

*’*  (5) 

"'Vm  '  V  “  ” 

These  two  assumptions  essentially  say  that  the  increment  in  velocity  at  any 
point  on  the  body  due  to  body  rotation  is  small  compared  to  the  translational 
velocity  of  the  body.  A  quick  check  on  thla  can  be  made  by  considering  a 
five  foot  long  body  traveling  at  sonic  speed,  say  1,000  fast  per  second.  If 
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the  center  of  rotation  Is  midway  along  the  body,  the  rotational  speed  necessary 
for  the  maximum  rotational  velocity  increment  to  reach  +10%  of  the  flight  speed 
is  given  by 


q  ■  0.1 


(xa  _  “  *.) 

1,(1)  • 


2.5  ft 


•  40  radians  par  second 
-  382  RPM  * 

Clearly,  it  ia  reasonable  to  aaaume  that  actual  rotational  speeds  will  not 
exceed  this  value.  An  additional  advantage  is  that  assumption  (5)  allows  the 
Integration  of  equations  (1-4)  to  be  carried  out. 

Carrying  out  the  integrations  of  equations  (1-4)  and  breaking  out  the 
linear  force  and  moment  terms  so  that  empirical  data  can  bs  used  Instead  of 
slender  body  values  gives  the  following  expressions. 

r,  •  -  7  \  '  ’"’-“I  '  *.,»>  +  V.,„] 

-  c  nc,  |  v  P  ♦  «2  •  I,  trv  ~ 

M  1  l  /J7T 

-  a  P  *  »21-  I.  -  .  (( 

J  /T7T  » 


F  -  -  7  P^MvaJ  CN  -  up  ur  IVtt  -  xa  m)  +  a2x.  1 
y  2  ••  o  N  ®  L  0  s,m  1  s,mj 

-  o  nc .  |  v  / v 2  4*  w2  •  I1  +  v 

'  1  [  tfTj- 


+  r  w2  •  I.  +  r  (rv  -  W)  .  i 

2  /2  .  2 "  3 

J  /v  +  w 


1-317 


.  l»*  M  t’ 

Ik  ?• 


M  -  \  PtB|u|waJ(xitin  -  x)^  -  np^uq  [e2tf  -  x^)2 

-  *ft« +  **]  ♦  {*  &***  ■  h 

♦  [» Xa •■■B’).  -  ,  42  +  ,2l.  i  - ,  .  ,  l 

I  JJ7T  J  JJ7T  > 


where 


N  "  "  f  “  X)CN  "  [*lU  ‘  x.,«)2 

-  •?*!,»  *  **]  *  |V  w2  '  !2 
J  irv+j*l+  r  /,2  +„2].  ,  i  | 

L  JJTT  J  5  JT7T  ’» 


<c 

ll  •  J  •  dx. 


4. 

22  •  j  ,(x..»  -  K.>dI'. 


■c  * 

h  "  I  *(x«,m  ’  x.)2dx*  *4  •  J  *2(x.,b  *  x.)dx. 

O  0 

l 

*5  *  |  a(xa,m  ’  x.)3dxa 


•  ■  radius  of  tha  cylindrical  portion  of  the  body 
0 

!0,  if  nose  ia  pointed 

a  ,  if  noaa  la  flat  facad 
o 

whara  n,  ca  ,  Cjj  and  x  ara  to  ba  determined  from  empirical  correlations.  To 
avoid  unnacSaaary  computar  run  time,  available  methods  ware  reduced  to  the 
aiinpleat  poasible  forma  which  did  not  sacrifice  accuracy  substantially3. 
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EQUIVALENT  ANGLE  OF  ATTACK  FORMULATION  FOR  FINS 


The  equivalent  angle-of-attaek  concept  la  daaorlbad  in  detail  in  References 
6  and  7.  In  brief,  the  idea  ia  to  calculate  aomahov  an  equivalent  angle  of 
attack,  a  ,  eo  that 


(11) 


where 

CN  ■  normal-force  coefficient  acting  on  fin  i  baaad  on 
i  planfora  area 

C.  ■  normal-force  coefficient  acting  wing  a^one  compoaad 
w  of  two  opposing  fine  with  aama  expn  d  plenform  as  fin  1 
joined  at  the  root  chorde.  Reference  area  la  the  plan- 
form  area. 


Equation  (11)  repreaanta  an  attempt  to  correlate  the  variations  of  the  fin 
normal  force  due  to  body  angle  of  attack,  body  roll  angle,  and  fin  deflection 
through  a  single  parameter,  a^n.  The  quantity  is  the  wing-alone  normal- 
force  coefficient.  The  wing  alone  ia  obtained  by  removing  the  body  between  two 
opposing  fine  and  joining  them  together  at  their  root  chords.  Uaing  experi¬ 
mental  values  for  Cnw  allows  the  incorporation  of  nonlinear  ef facta. 

Consider  the  side-edge  view  of  a  fin  shown  in  Figure  3.  The  first  step 
is  to  compute  the  velocity  components  Vpi  end  Vni  seen  by  the  area  centroid 
of  the  fin.  The  velocity  Increments  due  to  body  rotation  are  included. 
Component  Vpi  la  parallel  to  the  body  axis.  Component  Vni  la  perpendicular 
to  the  plane  the  fin  would  occupy  if  5*  were  aero. 

Due  to  flow  around  the  body  in  the  crossflow  plane,  the  average  normal 
velocity  aeen  by  pointa  on  the  fin  ic  lncreaaed.  This  phenomenon  la  known  aa 
Beskin  upwash4  and  la  primarily  a  function  of  the  shape,  angle  of  attack  and 
Mach  number  of  the  body7. 


Uaing  slender  body  theory,  one  can  show  that  the  ratio  of  the  normal  force 
acting  on  two  opposing  fine  in  the  presence  of  the  body  with  no  aidasllp  and 
cero  fin  deflection  to  the  normal  force  acting  on  a  wing  alone  composed  of  the 
two  fine  la  a  function  of  the  ratio  of  the  body  diameter  to  the  fin  open  only. 
This  ratio  la  called  Ky. 


For  the  present  work,  we  ehall  assume  that  the  effective  normal  velocity 
aeen  by  the  area  centroid  of  fin  i  ia  given  by  Ku'  when  the  fin  la  undeflected. 


Aa  a  firet  step  in  the  calculation  of  the  effects  of  fin  deflection,  we 
assume  that  fin  1  is  the  only  fin  on  the  body  (no  fin-fin  interference).  Then 
the  velocity  component  parallel  to  the  fin  root  ohord  is  given  by 


V’ 

pi 


V"  5i  -  4i 


(12) 
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The  velocity  component  perpendicular  to  the  fin  la  given  by 


v»i '  “''V" s* +  V‘" Sl 


To  include  vertex  effect.,  the  Increment  In  fin  normal  velocity  induced  by  the 
trailing  vorticea  from  the  wing  aection  ia  added  to  to  get 

v;  •  Vc/”*  si +  vp1[,ln  *i +  ,1"<‘“«fl>vi  ]  UA 

A  method  for  eatimating  (Ao^)^  l*  *lv*n  Ut#f  10  th*  PBper“ 

To  account  for  fin-fin  interference,  a  factor  A  ia  applied  to  the  deflec¬ 
tion  angle  of  aach  fin  aa  followa 

V’  ■  K..V  coa  5.  +  V  ain(4.A.)  +  *in(Aa  ) 
nt  1  piL  1 

-  I  ain(d  A  )  (l 

J-l  *  *  J 

SSSSSSSSriSi -Sa* 

alendar  body  theory  are  given  in  Reference  3. 

Tha  reeulta  of  Equationa  (12)  and  (15)  give  ®u  f°llOWB 


Then,  the  force  on  the  fin  ie  given  by 


where 


\  •  vv  'haW 


q  „  i  0  (v2  +  v2  ) 

4i  2P»  Pe  »4 


Methodology  for  computing  the  wing-alone  normal  force  coefficient  ia  given  In 
the  next  aection. 
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FIN  NORMAL  FORCE 


The  wing-alone  normal-force  coefficient  for  low  aspect  ratio  fins  has  been 
correlated  In  a  manner  similar  to  that  used  by  Allen  for  slender  bodies'4.  The 
equation  is 

"4  C.,  sin  2a  cos  a  +  K  sin^a,  0  x.  a  4  90®  (19) 

c  N 

w  a 

where  is  the  linear  normal-force  coefficient  slope  and  K  is  the  dreg  coef¬ 
ficient  8f  the  wing  when  it  is  normal  to  the  flow.  In  general  the  factor  K 
depends  on  aspect  ratio  and  Mach  number.  In  Reference  3,  correlations  for 
and  K  are  presented.  Also  presented  are  similar  correlations  for  high  aspect01 
ratio  swept  planforma  and  correlations  for  longitudinal  and  lateral  positions 
of  the  fin  center  of  pressure. 


WING-TAIL  INTERFERENCE 

The  method  used  here  for  the  computation  of  wing-tail  interference  is 
essentially  that  of  Reference  2.  However,  since  thut  methodology  was  derived 
for  an  unrolled  cruciform  or  monoplane  missile,  some  modification  is  needed 
to  handle  the  more  general  cases  encountered  for  this  work.  The  approach  here 
la  to  treat  the  missile  as  if  it  were  unrolled  but  with  the  crossflow  velocity 
equal  to  i/v'2  +  w-  and  with  the  wing  aection  developing  the  actual  normal  force 
component  which  la  parallel  to  and  in  the  seme  direction  aa  the  croaaflow 
velocity  vector  as  ssen  by  the  body.  The  methods  of  Reference  2  are  used  to 
compute  an  equivalent  angle  of  attack  acting  parallel  to  the  crossflow  velocity 
vector  at  the  axial  location  of  the  tail  fin  area  centroid.  That  angle  of 
attack  is  than  resolved  into  components  normal  tu  each  of  the  tail  fine  at 
their  actual  orientations.  The  msthematical  details  are  given  in  Reference  3. 


EXAMPLE  CASES 


SURFACE-TO-AIR  MISSILE 

For  the  surface-to-air  mlealle  simulation,  it  was  assumed  that  the  rocket 
motor  had  burned  out  and  that  the  guidance,  ordnance  and  autopilot/battery 
sections  had  been  separated  from  the  missile  as  shown  In  Figure  4.  The  piece 
of  debris  considered  was  tha  aft  fragment.  The  computed  static  margin  for  the 
fragment  was  found  to  be  negative  for  the  transonic  and  supersonic  speed  range. 
Hence,  the  fragment  can  be  expected  to  tumble  If  the  control  surfaces  are  not 
de  fleeted. 

Results  for  the  computed  trajectory  are  given  in  Figure  4,  Initially,  the 
miaaile  Is  In  a  10°  dive.  At  t  ■  t„,  the  forebody  is  separated.  The  aft 
fragment  is  assumed  to  be  given  a  slight  initial  q  by  tha  separation  event. 

It  immediately  pitches  up  and  rapidly  tumbles.  The  velocity  of  the  fragment 
drops  quickly  and  the  aft  fragment  hita  the  ground  roughly  2,000  feet  short  of 
the  aiming  point.  It  has  been  assumed,  of  course,  that  the  fragment  does  not 
disintegrate  during  the  hlgh-q  pitch  up.  It  is  interesting  to  note  thut  the 
trajectory  does  not  appear  to  be  ballistic  until  t  -  t  »  2  seconds. 
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Essentially,  this  means  that  downrange  travel  of  the  fragment  would  probably 
have  been  underestimated  by  n  two-dimensional  computation  based  on  an  average 
drag  coefficient. 


RYAN  BQM-34A  TARGET 

The  geomatrical  characteristics  of  the  BQM-34A  target  as  given  in  Refer¬ 
ence  8  are  shown  in  Figure  S.  The  debris  fragment  to  be  considered  is  that 
portion  of  the  configuration  which  is  aft  of  tha  dashed  line  in  the  side  view 
(nose  and  engine  gone).  Because  of  the  methodology  limitations  it  was 
necessary  to  model  the  body  aft  of  the  fragment  as  a  circular  cylinder.  Three 
trajectories  were  run  for  the  BQM  aft  fragment  to  illustrate  the  effect  of 
small  changes  in  the  predicted  static  margin.  At  the  time  of  separation,  tha 
vehicle  la  assumed  to  be  in  level  flight  at  M,  ■  0.75  at  an  altitude  of  500 
feet  AGL.  The  computed  results  era  shown  in  Figure  6.  During  trajectory  it  1, 
the  damaged  vehicle  pitched  up  to  about  1.7*  angle  of  attack  and  maintained 
that  attitude  until  Impact.  During  trajectory  it 2,  the  vehicle  quickly  pitched 
up  to  1.5*  angle  of  attack.  It  continued  to  slowly  Increase  the  pitch 
attitude  until  the  top  of  its • trajectory  at  which  point  the  angle  of  attack 
was  approximately  40*.  It  then  nosed  over  and  fail  to  the  ground.  During 
trajectory  (73,  the  vehicle  slowly  pitched  up  to  roughly  Is  angle  of  attack. 

At  t  -  t0  *  6  seconds  (4,000  feet  downrange),  it  no  longer  had  sufficient  speed 
to  maintain  level  flight  and  descended  until  Impact. 

Thu  "trimmed"  flight  behavior  of  the  BQM  aft  fragment  as  shown  in  Figure  6 
appears  to  be  a  result  of  the  tendency  of  the  center  of  pressure  of  the  frag¬ 
ment  to  move  aft  as  tha  angle  of  attack  is  increased.  Hence,  it  is  possible 
for  there  to  be  an  angle  of  attack,  a0,  such  that  for  a  <  a0  the  fragment  is 
unstable.  Any  alight  disturbance  would  cause  the  vehicle  to  pitch  to  a0.  This 
phenomenon  causes  the  particular  flight  behavior  encountered  to  be  very  sensi¬ 
tive  to  the  center-of-presaure  and  center-of-gravlty  locations  if  the  stability 
of  the  fragment  at  small  angle  of  attack  is  nearly  neutral. 


SUMMARY 

A  set  of  methods  has  been  developed  for  computing  the  longitudinal,  lateral 
and  control  aerodynamic  characteristics  of  a  wide  range  of  missile  debris 
fragments.  The  methods  have  been  Incorporated  into  a  computer  program  which 
simulates  the  six-degree-of- freedom  trajectories  of  the  fragments.  Sample 
casea  presented  in  this  paper  and  in  Reference  3  demonstrate  that  ballistic, 
tumbling,  coning  and  "trimmed"  flight  trajectories  cen  sll  be  predicted  for 
reasonebie  cost.  Typical  running  times  range  roughly  from  0.5  to  5  times  real 
time  on  a  CYBER  175  computer  depending  on  the  complexity  of  the  motion.  While 
the  computer  program  waa  designed  to  determine  the  lethality  of  a  missile  once 
it  has  broken  up,  it  can  also  be  usad  for  range  safety  studies. 
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LIST  OF  SYMBOLS 


a 


local  radius  of  body,  feet 


a 

o 

al 

cNi 

CNW 


radius  of  cylindrical  portion  of  body  alona,  fast 
radius  of  body  at  xf  -  0,  fast 

coafficlant  of  normal  fores  acting  on  fin  i{  normal  force/QaSr#£ 

normal-fores  coafficlant  of  wing  alone  formed  by  putting  together 
two  opposing  fins  at  their  root  chords;  normal  force/QaSra£ 

derivative  of  normal-force  coefficient  with  respect  to  angle  of 
attack,  a,  at  a  ■  0 


c,jc  crossflow  drag  coefficient 

FNi  normal  force  acting  on  fin  1,  lbf 

F  ,F  ,F  component  in  body-fixed  coordinate  system  of  force  acting  on 
y  f  ragmen tr  lbf 

K  drag  coefficient  of  a  wing  when  it  is  normal  to  flow 

1C.  ratio  of  normal  force  acting  on  two  opposing  fins  In  presence  of 

body  to  normal  force  of  wing  alone  at  same  angle  of  attack  as  body; 
no  sideslip  and  no  fin  deflection 


L,M,N 


l 


components  In  body-fixed  coordinate  system  of  moment  acting  on 
fragment,  ft*lbf 

length  of  body,  feet 


P.q.r 

M.t 


u,v,w 


u,v,w 


components  along  body-fixed  coordinates  of  rate  of  rotation  of  body 
about  its  center  of  mass,  radian/sec 

2 

rate  of  change  of  p,  q,  r  with  respect  to  time,  radian/sec 

1  p-(Vo  *  Vn  > 

2  Pi  ni 

components  along  body-fixed  coordinates  of  velocity  of  body  center 
of  mass,  ft/sac 

o 

rate  of  change  of  u,  v,  w  with  respect  to  time,  ft/sec 


Pi 


velocity  component  parellel  to  body  axis  at  fin  i  area  centroid 
with  no  fin  deflection,  ft/eec 


velocity  component  normal  to  fin  i  at  fin  i  area  centroid  with  no 
fin  deflection,  ft/sec 
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X  I 


LIST  OF  SYMBOLS  (Concluded) 


x 

s 

X 

s,m 


P 


DO 


rearward  distance  along  body  axis  from  nose  tip,  ft 

distance  from  nose  tip  of  fragment  to  center  of  mass,  ft 

axial  location  of  center  of  pressure  measured  from  nose  tip  of 
fragment,  ft 

angle  of  Incidence,  t«n~*(VC];/u) 
equivalent  angle  of  attack  of  fin  1 
deflection  angle  of  fin  i 

parameter  accounting  for  finite  length  of  body,  dimensionless 

fin-fin  interference  factor  due  to  deflection  of  fin  i 

increment  in  equivalent  angle  of  attack  of  tail  fin  i  due  to 
presence  of  vortices  shed  from  wing  asction 

missile  roll  angle,  tnn-1(w/v) 

3 

atmospheric  density,  slugs/ft 
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(a)  Nose  blown  offi  -  a.l  Nosaa  of  varying  slsndernesa 

a. 2  Cylindrical  bodlea  with  or  without 
winga  and  talla 


(b)  Tall  blown  offt  -  b.l  Noaaa  with  aftarbodlaa  of  varying 

alenderneaa,  with  wing  aaction 


b.2  Short  cylindrical  aaction  with 
tall  aaction 

a 

© 

a 

(c)  Noae  and  tail  blown  off:  -  c.l  Noses  of  varying  slenderness 

c.2  Cylindrical  sections  of  varying 
length,  with  wing  section 

c.3  Cylindrical  bodies  with  wing 
and  tail  sections 


Figure  1.  Possible  debris  fragment  shapes 
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Figure  2.  Coordinate  system  fixed  in  debris  fragment 
and  used  in  force  and  moment  calculation 


Plane  of  fin  i 


Leading  edge 
of  fin 


Line  parallel 
to  body  axis 


.Trailing  edge 
S  of  fin 

'  \ 
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Figure  3.  Side-edge  view  of  fin  i  showing  velocity  components 
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